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Abstract 

We present exact explicit analytical results describing the exact ground state of four electrons 
in a two dimensional square Hubbard cluster containing 16 sites taken with periodic boundary 
^ ', conditions. The presented procedure, which works for arbitrary even particle number and lattice 
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sites, is based on explicitly given symmetry adapted base vectors constructed in r space. The 
Hamiltonian acting on these states generates a closed system of 85 linear equations providing by 



'^ . its minimum eigenvalue the exact ground state of the system. The presented results, described 

O I with the aim to generate further creative developments, not only show how the ground state can be 



exactly obtained and what kind of contributions enter in its construction, but emphasize further 



^ . characteristics of the spectrum. On this line i) possible explications are found regarding why weak 

coupling expansions often provide a good approximation for the Hubbard model at intermediate 

cn 

^r I couplings, or ii) explicitly given low lying energy states of the kinetic energy, avoiding double 



occupancy, suggest new roots for pairing mechanism attracting decrease in the kinetic energy, as 
emphasized by kinetic energy driven superconductivity theories. 

PACS numbers: 



I. INTRODUCTION 



The growing interest in systems of highly correlated electrons has proposed many theoret- 
i cal descriptions where the dominant effects are treated in the frame of the Hubbard model 



Hubbarc 



Lieb and Wu 



1963 ). Since the exact solution of the model is known only in one dimension 



1968^ . and the majority of driving forces in the field, as superconductivity 



in cuprates, are two dimensional (2D 



have been proposed (jBaeriswvl et al. 



prob lems, a great variety of approximate treatments 



19951 ) in order to accommodate a suitable theoreti- 



cal framework. Despite several years of intensive studies, it is apparent that the nece ssary 



theoretical skills and tools to deal with this problem are still in fact relatively poor (JLieb 
19941) . 



While most of the ongoing effort relating the Hubbard model is being concentrated 
on the general case of large electron density, more than ten years ago it has been recog- 
nized that the less an alysed low-density limit not only that retain main aspects related to 



the model behaviour (^Falicov and Proetto 



1223, 



Papavassiliou and Yartse 






19921 ). but can 



provide key knowledge which cou ld drive further at least nonperturbative developments 



( Fabrizio. Parola and Tosatti 



19911 ). firstly because is more hope to find non- approximated 



solutions holding by their nature essential information. The history of this background has 



started with the exact solution of the two-e 



torus (IChen and Mei 



( Parola et al. 



m 



198G 



Mei and Chen 



ectron ground state problem on an arbitrary large 



1988^ . solution of the two electron problem in 2D 



199flD . and continu ed with the three electron problem extended nonpertur 



batively in the low density limit (JFabrizio. Parola and Tosattilll99l[ ). The numerical treat- 
ment of the four electron problem on four sites follows by the test of the Bardeen-Cooper- 
Schrieffer and re sonant- valence-bond wave functions as approximated ground states of the 
Hubbard model (Falicov and Proettdll993l ). test of the unrestricted Hartre e-Fock solution 
(ILouis et al.lll992l) . the calculation of Huckel-Hubbard correlatio n diagrams (IZhu and Jiane 
1993[) . and deduction of absorption spectra of TCNQ particles flPapavassiliou and Yartsev 
992|). The concurrently made numerical d evelopments on 4 x 4 clusters around half filling 



( Galan and Vergei 



1991 



Parola et al. 



19911) corrobor ated by group-theoretical studies of the 
4x4 square system (iFano. Ortolani and Parolalll992l l also lead to the de velopment of the nu- 



merical des cription of the 4 particle problem for energy level statistics (|Bruus and d'Auriac 



mm, 



19971 ). In the last years the study of symmetry properties for Hubbard clusters in 



holding 4 electrons has been given (jAcquarone et al. 



general has continued ()Tieniberdll998^. extensions to Hubb ard-Holstein model for clusters 



m 



19990, ^^^ fo^ spinless fermion case 



num erical characterization of the 4 particle problem goes up to clusters of 20 x 20 exten- 



sion (jZhang and Henlev 



2004f ). Despite the invested efforts, exact analytical results holding 



essential information are not known at the moment in the problem. 

The last period has provided strong new motivations for the further development of 
the study of the low concentration limit of the Hubbard model. The reason for this 
is that experimentally one started to encounter in condensed matter context rapidly in- 
creasing situations containing small numbe r of particles confine d in a system or device, as 
for example in the case of quantum dots (IMaksvm et £11.1120001) . quantum well structures 



Kochereshko 



et al. 



Sackett et al. 



] 



20031 ) ■ mesoscopic systems (|Halfpapl 120011 ). experimental entanglement 



20001 ) etc. Between the studied experiments, seve ral are directly co nnected to 



the Hubbard model, as in the case o 



dots (JBusser. Moreno and Dagotto 



charge transfer complexes (ISaito et al. 



20041 ). or mesoscopic grains (IBovaci. Gedik and Kulik 



20011). quantum 



20001 ). Furthermore, it has been observed that several measurements on strongly correlated 
systems are acceptable reproduced on small H ubbard cluster s, as in the case of the X-ray ab- 
sorption for nickelates studied on 4-site cluster ( Okada 12004 ). The same can be stated for the 



4 X 4 cluster case used in describi ng manganites (IHelberd 



20021 ). photoemission intensities (JEroles. Batista and Aligia 



2001 



Srinitiwarawong and Gehring 



19991), or thermodynamic prop 



erties provided by processes in the vicinity of the Fermi surface (jChiappe et al 



mM) 



As can be observed, the accumulated knowledge to the present date regarding the ex- 
treme low density limit described by few particles present in the system, excepting the 
one and two particle cases, means only symmetry properties, and numerical results (often 
approximated), which usually conceal explicit properties which could generate creative ad- 
vancement at the level of the theoretical description and deep understanding. This is a 
regrettable and unfortunate situatio n, since as is known from the accumulated experience in 



solving e xactlv manv-bodv svstems (|Mattis 



stressed (IFabrizio. Parola and Tosatti 



imi 



1993 K or a s several times has been accentuately 
lLieblll994[ ). key aspects of the unapproximated 



descriptions are often hidden in the few particle cases. 

On the presented background, with the aim to fill up at least partially this gap for 
the two dimensional Hubbard case, and being driven by the intention to provide essential 
information for the subjects mentioned above, we present the exact ground state of four 



interacting electrons in a 4 x 4 Hubbard cluster. Our main purpose is not to hide potential 
essential characteristics behind numerical results, but to present explicit expressions, basic 
resulting properties, and visible characteristics, which we strongly hope, are able to polarize 
the creative thinking, and advancement in the field. 

Our results are based on an r space description, which in our opinion is unjustly over- 
shaded in the last period in treating such problems, in contradiction with its deep ability 
to bring to light essential characteristics. We construct first a base vector set starting from 
symmetry properties, and then show how a closed system of linear equations containing 85 
components can be constructed characterizing the ground state manifold and providing by 
its secular equation, through its minimum eigenvalue, the ground state wave function and 
ground state energy. After this step the properties of the ground state are analyzed. The 
concretely described situation is an L x L = A^ = 16 square lattice with periodic bound- 
ary conditions in both directions, and Np = 4 particles. But the method itself works for 
arbitrary even number of particles and arbitrary even L. The explicitly presented base vec- 
tors and equations providing the ground state, not only show how the ground state can be 
constructed and what kind of components build it up, but provides an insight into other 
properties of the spectrum as well. Especially the characteristics regarding the kinetic energy 
eigenstates should be mentioned on this line, important for understanding aspects related 
to perturbative expansions, or kinetic energy driven superconductivity. 

The remaining part of the paper is structurated as follows. Section II. presents the Hamil- 
tonian, the deduction procedure and the ground state wave functions. Section III. describes 
physical properties of the deduced eigenstates. Sect. IV. presents the conclusions of the 
paper, while the Appendices A - B presenting mathematical details, close the presentation. 

II. HAMILTONIAN AND GROUND STATE WAVE FUNCTIONS. 
A. Presentation of the Hamiltonian 

The Hamiltonian we use has the form of a standard Hubbard Hamiltonian 

<i,j>,a i 

where c]^ creates an electron at site i with spin a, t is the nearest-neighbour hopping 
amplitude, U is the on-site Coulomb repulsion, and < i,j > represents nearest-neighbour 



sites, taken into account in the sum over sites only once. During our study we consider 
an L X L cluster, N = L"^ , with periodic boundary conditions in both directions. All the 
presented results will relate H /t, the only one microscopic parameter of the problem being 
u = U/t. The case A^ = 16 will be presented in details, although the method is applicable 
for arbitrary even A^. Since Ap = 4, the band filling for the presented case is 0.125. We 
further mention that we consider during this paper u > 0. 

B. The construction of the base wave vectors 

1. The basic wave vector elements 

In order to describe the Hilbert space region containing the ground state for the Ap = 
4 particle problem, we use an r-space representation for the wave vectors. In order to 
characterize this, first the lattice sites of the considered system are numbered starting from 
the down-left corner, the numbering being given along the lower lattice sites line, then going 
upward with the notation, as shown in Fig.l.a. In the considered system (taking into account 
that the ground state is a singlet state), three type of particle configurations may occur as 
depicted in Fig.l.b,c,d. 

1) We can have two double occupancies at sites i and j represented by two dots at sites 
i and j (see Fig.l.b). 2) We may have a double occupancy at site i and two electrons with 
opposite spins at sites j and k. In this case, the double occupancy at site i it is denoted as 
before by a dot at i, while the two electrons with opposite spins placed on the different sites 
{j, and k) are depicted by a dotted line connecting the sites j and k (see Fig.l.c). Finally, we 
may have four single occupancies placed on different sites. From these, the electrons placed 
at i and j have spin a, while the electrons at sites k and / have spin —a. This situation will 
be represented by two continuous lines connecting the sites i and j, and k and /, respectively 
(see Fig.l.d). The notations i,j, k, I are representing the numbering of the sites as specified 
in Fig.l.a. 

The mathematical expressions connected to the states represented in Fig.l.b,c,d are given 
as follows. The two double occupancies present in Fig.l.b provide the state 

\b) = (iAMAi)\^)^ (2) 

where |0) represents the bare vacuum with no fermions present. The one double occupancy 



and the two electrons with opposite spins depicted in Fig.l.c are described as 

\c) = iiAimiAi) + (i,Ai)m ■ (3) 

For the here presented situation we must have j ^ k and i ^ j,i ^ k respectively. Finally, 
in describing mathematically the four single occupancies presented in Fig.l.d, we have 

\d) = \A,AMAi) + AiiMA,im^ (4) 

where for the mathematical clarity, j > i and I > k is required (all sites being considered 
different). 

2. The translation of the basic elements 

In order to construct the base vectors, the basic wave vector elements are translated by 
the operation T. The T operator is considered a linear operator, hence the relation 

T{A + B)=T{A)+T{B), (5) 

holds, where A and B represent particle configurations as depicted in Fig.l.b,c,d. Further- 
more T{A) means that the configuration of particles represented by A is translated to each 
site of the system (for A^ = 16, the translation is made 16 times), and the obtained contribu- 
tions are all added (see Fig. 2). The argument of the T operator contains usually four terms 
since the starting particle configuration is rotated by 180 degrees along the x, y, and z axes, 
and then the obtained contributions are added. Less than four terms in the T argument 
means that the mentioned rotations provide the same starting particle microconfiguration. 
The translations are such effectuated that the relative inter-particle positions are all main- 
tained, and arriving near the border of the system, the presence of the periodic boundary 
conditions are explicitly taken into account. The translations must be given at the level of 
the graphic plots, all the obtained 16 graphics must be added (see Fig. 2), then each graphic 
must be transformed in a mathematical form according to the rules presented in Eqs.(j2]IH). 
At this point we must underline, that when the mathematical expression is written for a 
given graph, a negative sign may emerge in some cases given by the restrictions prescribed 
for Eqs.(j2]|^. For example, the last two terms in both rows of Fig.3 have mathematical ex- 
pressions with negative sign. Similar reasons lead to sign changes obtained after rotations. 



reason which generates the sign changes observed in some T(...) arguments of base vectors 
presented in Figs. 4-18, for example in the case of the second contributions of |10), or |11), 
etc. Furthermore, in some cases (see for example |12) or |15)), the rotations provide the 
starting microconfiguration of particles, which leads to the decrease of the number of terms 
present in the argument of the T operation. The sign changes also lead in some cases to 
complete elimination of some possible particle microconfigurations. 

3. The base wave vectors 

Starting from the rules and operations described in details in Sec. II. B. 1-2, the 85 base 
wave vectors can be explicitly constructed as given in Fig.4-18. These are denoted by \n), 
n = 1,2, ...,85, and are all orthogonal vectors. As can be seen, the base wave vectors are 
obtained by subtracting two type of contributions. The first part of these is obtained from 
a starting particle microconfiguration (presented as first terms in the argument of the first 
operator T) which is rotated by 180 degrees along the x, y and z axes, and translated by the 
operation T after this step. The second (subtracted) part of the wave vectors is obtained 
from the first part by a rotation of the elements by 90 degree along the z axis. 

C. The ground state wave function 

First of all we are constructing a closed system of equations containing the base vectors 
described above. 

This system of 85 equations is obtained as follows, i) We start from the most condensed 
and most interacting particle microconfiguration depicted by the base vector |1) (see Fig. 4), 
which cannot be eliminated on physical grounds from the ground state, ii) By applying the 
Hamiltonian on the base vector |1) (see the first equation of Appendix A), we obtain as a 
result new base vectors (|2) and |3)), which have the same symmetry properties as |1). iii) 
Since the Hamiltonian does not change the symmetry properties of the base vectors, one 
must further apply H on the each resulting new base vector, up to the moment in which 
the obtained system of equations closes. This happens for the studied system at the 85th 
equation. The obtained closed system of equations is explicitly presented in Appendix A. 
The minimum energy solution of the presented system of equations represents the ground 



state of the system. 



III. PROPERTIES OF THE SOLUTIONS 



A. The check of the ground state energies 

First of all one can check that the system of equations presented in Appendix A indeed 
provides the ground state of the problem. The obtained minimum eigenvalues provided by 
Appendix A are 



u = 0.0, E 

u = 0.5, E 

M = 1.0, E 

M = 1.5, E 

M = 2.0, E 

M = 2.5, E 

M = 3.0, E 

M = 3.5, E 

M = 4.0, E 



-12.0000000000000000 
-11.9126285145094126 
-11.8364690235642218 
-11.7695877912829268 
-11.7104580743242632 
-11.6578605150913841 
-11.6108110503797608 
-11.5685079573602838 
-11.5302924026297138 



(6) 



which perfectly coincide with the numerically exact minimum eigenvalue s obtained from ex- 



act n umerical diagonalization in the 14 400 dimensional full Hilbert space (|Sorella and Hubsh 
20001 ). We note that the eigenvalues are given in t units. 



B. The U independent eigenstates 



The study of the solutions of the system of equations presented in Eq. ()Al|) leads to an 
extremely interesting observation, namely that almost half (exactly 40%) of the eigenstates 
are U independent, consequently are eigenstates of the non- interacting system as well. For 
example, the U independent 18 eigenvectors corresponding to zero eigenvalue are presented 
in Appendix B. Besides these, one has also 16 eigenf unctions corresponding to eigenvalues 
±4, ±8, which are also U independent (see exemplification in Appendix B). The presence of 
such eigenstates is important for following reasons. 



8 



First of all, several authors observed (JGalan and Verged ll99lL JMetzner and VoUhardt 



19891) that weak coupling expansions often provide a good approximation for the Hubbard 
model at intermediate coupling. A possible explication for this is the emergence of a huge 
number of eigenstates in the spectrum of the Hubbard model with non-zero interaction, 
which are present in the non-interacting case as well in the spectrum of the model. 

A second aspect which must be mentioned here is that the eigenstates related to zero 
eigenvalue presented in Appendix B are in fact many-body eigenstates of the kinetic energy 
term with no double occupancy. Such states apparently totally avoid double occupancy at no 
cost of ener gy, consequentlv are important in the studv of the pairing mechanism in the Hub 
bard model fjCini. Perfetto and Stefanuccill2001 , 



Cini and Stefanucciii2001 



Perfetto and Cini 



2004^ . and as seen from the presented results, easily emerge in the spectrum. 

At this step we underline that contr ary to the practice used in the theoretical studies up 



to this moment (Balzarotti et al. 



2004^ . double occupancy avoiding eigenstates emerge not 



only at zero energy, but also at much lower energy values, closely situated to the ground 
state energy. We exemplify this statement by the eigenstate \v2) presented in Eq. flB2j) . 
corresponding to energy —8 (in t units). These eigenstates being U independent, their 
energy remain unchanged if the strength of the on-site interaction is increased, while the 
ground-state energy increases if U is increased. Consequently, starting from these states, 
the contribution in the pairing mechanism of the kinetic energy eigenstates not containing 
double occupancy could be much more efficiently taken into account than considered up 
today. We note that this problem has interconnections to the problem of the kinetic en- 
ergy driven superconductivity as well, strongly debated in the literature in the last period 



(Anderson 



1995 



Feng 



2003 



irschi 



Hirsc 



2004, 



observed (JEckl. Hanke and Arrigoni 



2003 



Yokovama et al. 



200- 



van der Marel et al 



). Ind eed, experimentally is 



20031) that especially in Bi 



based cuprates, a decrease in the kinetic energy is observed at the su perconducting transi- 
tion, which contradicts the usual (for example BCS) pairing theories ( Hirschll2004r i. If the 
pairing process could be described on the states of the type |f2), exemplified in our knowl- 
edge in exact terms for the first time in this paper, a such a decrease could be much better 
understood. 



C. The U dependent eigenstates 

The major part of the spectrum (60%) is U dependent. From the eigenstates entering in 
this category, there are states whose eigenvalue is exactly u, for example |i(;) = |8) — |9) + 
|39) + |40) — |41). Such states are interesting for two reasons. First, a detection possibility of 
such states would lead to a direct experimental measurement of the U/t ratio. Second, since 
all the contributions in \w) have double occupancy d = 1, this vector is also an eigenstate 
of the interaction term with eigenvalue u, and in the same time, the eigenvector of the 
kinetic energy term with eigenvalue zero. Consequently, the study of the zero kinetic energy 
eigenvalue states must be given with care, since such type of states, as exemplified before, 
not necessarily are avoiding the double occupancy. 

The remaining part of the deduced U dependent states, are not eigenstates of the kinetic 
energy, and the interacting part of H separately, but are eigenvectors only for the sum 
of both, e.g. the whole Hamiltonian. As seen from Eq.®, the corresponding eigenvalues 
have a smooth (less than linear) U dependence, and the ground states is obtained from this 
category. 

D. The system generating the ground state 

In order to enhance further developments and enlighten further creative thinking, we pro- 
vide explicitly in Appendix A the equations generating the ground state, together with the 
explicit form of the base vectors providing the ground state (see Figs. 4-18), hence describing 
the Hilbert space region in which this is placed. We mention that if the base vectors are 
defined only through the operation T, without subtracting two T components as shown in 
Figs. 4-18, a system containing 176 equations arises. This however can be cast in two block 
diagonal components, one of which provides the equations presented in Appendix A. 

IV. SUMMARY AND CONCLUSIONS 

Driven by the aim to provide explicit expressions generating creative developments, we 
present exact analytical results describing the ground state of the two dimensional Hubbard 
model taken on an L x L = A^ square lattice at A^ = 16 (and periodic boundary conditions 
in both directions), containing Np = 4 particles. The presented procedure allows the de- 

10 



scription of an arbitrary even L and Np, and is based on symmetry adapted base vectors 
constructed in r space. The Hamiltonian acting on the described base vectors provides a 
closed system of hnear equations (whose number is 85 for A^ = 16 and Np = 4), leading by 
its secular equation, through its minimum eigenvalue, to the ground state wave function and 
ground state energy of the system, deduced in a Hilbert subspace with almost three orders 
of magnitude smaller in dimensions than the full Hilbert space of the problem. The deduced 
explicit eigenstates characterize also other properties of the spectrum: i) The large number 
of eigenstates which remain eigenstates of the non-interacting Hamiltonian as well shows 
why weak coupling expansions often provide a good approximation for the Hubbard model 
at intermediate coupling, ii) Zero energy eigenstates of the kinetic energy term {Hkin) which 
are eigenstates of the Hamiltonian as well show how energy increasing double occupancies 
can be avoided providing a possible support for the kinetic energy driven superconductivity, 
iii) Low energy eigenstates of the kinetic energy term which completely avoid double occu- 
pancy emphasize potentially new pairing possibilities in the low energy part of the spectrum 
in the context of the kinetic energy driven superconductivity, iv) Zero energy eigenstates 
of the kinetic energy term corresponding to double occupancy one underline that the zero 
energy H^in eigenstates must be handled with care, since as exemplified, these states not 
necessarily represent double occupancy avoiding states. 

Similarly obtained solutions for higher L or Np, corroborated with the study of the 
emerging changes in the system of equations describing the ground state, remain a challenge 
for future developments. 
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APPENDIX A: THE LINEAR SYSTEM OF EQUATIONS CONTAINING THE 
GROUND STATE. 

This Appendix presents the 85 hnear equations describing the ground state of the studied 
system. 

iJ|l) = 2M|l) + |2)-|3), 

H\2) = |1) + u\2) + 8|4) - |6) - |7) + |10) + |11) + 4|12), 

H\3) = -8|1) + m|3) - 8|5) + |6) + |7) - 2|8) - 2|9) - |10) - |11) - 2|13) - 2|14) - 4|15), 

^|4) = |2) + 2w|4)-|17), 

^|5) = -|3) + |23), 

H\6) = -2|2) + |3) + u\6) - |16) + 2|17) - 2|18) - |19) + |20) + 2|21) - |23) + |24) 

+ |25) + |26) + |29) + |31), 
H\7) = -2|2) + |3) + u\7) - |16) + 2|17) + |20) + 2|22) - |23) + |25) + |26) + |28) + |30), 
H\8) = -|3) + u\8) - |16) - |21) - |22) + |23) - |27) - |33) - |34), 
^|9) = -|3) + u\9) - |16) - |19) - |20) + |23) - |24) - |27) - |32), 
if|10) = 2|2) - |3) + |16) - 2|17) - 2|18) - |19) + |23) + |24) - |25) - |26) - |28) 

-|31)- |32)-2|33), 
H\ll) = 2|2) - |3) + |16) - 2|17) + |23) - |25) - |26) - |29) - |30) - |32) - 2|34), 
H\12) = |2)- |25), 

H\13) = -|3) - |16) + |19) + |23) + |24) - |27) - |28) - |29), 
H\U) = -|3) - |16) + |23) - |27) - |30) - |31), 
^|15) = -|3)-|26)-|27), 
H\16) = -|6) - |7) - 2|8) - 2|9) + |10) + |11) - 2|13) - 2|14) + m|16) - 8|36) 

+4|46)-4|54), 
if |17) = -8|4) + |6) + |7) - |10) - |11) + u\l7) + 4|36) - 2|40) - 2|41) - 4|44) 

-4|47) + 2|63) + 2|64), 
i;^ |18) = -|6) - |10) + u\18) + |42) - |56) + |58), 
if |19) = _|6) - 2|9) - |10) + 2|13) + m|19) - 2|35) - 4|39) - 2|42) - 2|55) + |57) + |58), 
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^|20) = |6) + |7) - 2|9) + m|20) - 2|35) - 2|38) - 4|40) - 2|42) - |43) - |50) - |51) 

+ |57) + 2|61) + 2|65), 
H\2l) = |6) - |8) + u\2l) - 2|41) + |42) - |43) - |50) - |56) - |59), 
^|22) = |7) - |8) + m|22) - 2|38) - 2|41) - |51) - |59), 
iJ|23) = 8|5) - |6) - |7) + 2|8) + 2|9) + |10) + |11) + 2|13) + 2|14) + 2|37) 

+2|38) + 2|43) + 8|44) - 4|45) - 4|46), 
H\24:) = |6) - 2|9) + |10) + 2|13) - 2|35) - 4|39) - 2|55) - 2|56) + |57) - |58), 
if |25) = |6) + |7) - |10) - 111) - 8|12) + 4|45) + 4|48) + 4|54) - 2|59) - 2|60) 

+2|61) + 2|62), 
H\26) = |6) + |7) - |10) - 111) - 4|15) - 8|47) + 4|48) - |49) - |50) - |51) - |52), 
if |27) = -2|8) - 2|9) - 2|13) - 2|14) - 4|15) + 4|45) - |49) - |50) - |51) - |52) 

-8|53)-4|54), 
^|28) = |7) - |10) - 2|13) + 2|35) - 2|38) + 2|42) - |43) - |51) - |52) - |57) 

+2|62)+4|63)-2|65), 
iJ|29) = |6) - 111) - 2|13) + 2|35) - 2|37) - |43) - |49) - |50) + 2|56) - |57) 

+2|62)+4|63)-2|65), 
^|30) = |7) - 111) - 2|14) - 2|37) - 2|38) - |49) - |51) - 2|60) + 4|64), 
H\31) = |6) - |10) - 2|14) - 2|43) - |50) - |52) - 2|60) + 4|64), 
^|32) = -2|9) - |10) - 111) - 2|35) - 2|37) - 4|40) - |43) - |49) - |52) - 2|56) 

+ |57) + 2|61) + 2|65), 
H\33) = -|8) - |10) - 2|41) - |42) - |43) - |52) + |56) - |59), 
H\34) = -|8) - 111) - 2|37) - 2|41) - |49) - |59), 
iJ|35) = -|19) - |20) - |24) + |28) + |29) - |32) + |74) - |75), 
H\36) = -|16) + |17) + 2m|36) + |66), 

iJ|37) = |23) - |29) - |30) - |32) - 2|34) + 2|67) - |70) + |71), 
^|38) = -|20) - 2|22) + |23) - |28) - |30) + 2|67) - |70) + |71), 
H\39) = -|19)- |24)+m|39), 
^|40) = -|17) - |20) - |32) + m|40) - |66) + |67) - |69), 
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H\Al) = -|17) - |21) - |22) - |33) - |34) + u\Al) - |66) + |67) - |69), 

H\A2) = |18) - |19) - |20) + |21) + |28) - |33) + m|42) + |75) - |76), 

H\A3) = -|20) - 2|21) + 2|23) - |28) - |29) - 2|31) - |32) - 2|33) + 4|67) - 2|70) + 2|71), 

H\U) = -|17) + |23) + |67), 

^|45) = -|23) + |25) + |27) + |70), 

H\4Q) = |16) - |23) - |71), 

H\47) = -|17)-|26)-|69), 

H\48) = |25) + |26)-|72), 

H\m) = -|26) - |27) - |29) - |30) - |32) - 2|34) - 2|68) - 2|69) - |70) + |72) - |73), 

^|50) = -|20) - 2|21) - |26) - |27) - |29) - |31) - 2|68) - 2|69) - |70) + |72) 

-|73) + |74)-|75)-2|76), 
^|51) = -|20) - 2|22) - |26) - |27) - |28) - |30) - 2|68) - 2|69) - |70) + |72) - |73), 
iJ|52) = -|26) - |27) - |28) - |31) - |32) - 2|33) - 2|68) - 2|69) - |70) + |72) 

-|73)-|74) + |75) + 2|76), 
H\53) = -|27)- |68), 
H\54) = -|16) + |25) - |27) - |73), 
^|55) = -|19)-|24) + |74) + |75), 

H\56) = -|18) - |21) - |24) + |29) - |32) + |33) + |74) + |76), 
H\57) = |19) + |20) + |24) - |28) - |29) + |32) - |74) - |75), 
H\58) = 2|18) + |19) - |24) - |74) + |75) + 2|76), 
iJ|59) = -|21) - |22) - |25) - |33) - |34) - |70) + |72) + |73), 
^|60) = -|25) - |30) - |31) - |70) + |72) + |73), 
iJ|61) = |20) + |25) + |32) + |70) - |72) - |73) - |74) - |75), 
H\62) = |25) + |28) + |29) + |70) - |72) - |73) + |74) + |75), 
^|63) = |17) + |28) + |29) + |66) - |67) + |69), 
H\64) = |17) + |30) + |31) + |66) - |67) + |69), 
H\65) = |20)- |28)- |29) + |32), 
^|66) = 4|36) - 2|40) - 2|41) + 2|63) + 2|64) + u\66) + 4|78) - 4|85), 
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iJ|67) = 2|37) + 2|38) + 2|40) + 2|41) + 2|43) + 4|44) - 2|63) - 2|64) + 16|77) 

-4|78)+4|79), 
iJ|68) = -|49) - |50) - |51) - |52) - 4|53) - 8|80) + 4|81), 
^|69) = -2|40) - 2|41) - 4|47) - |49) - |50) - |51) - |52) + 2|63) + 2|64) - 4|79) 

-8|80)-4|85), 
H\70) = -2|37) - 2|38) - 2|43) + 4|45) - |49) - |50) - |51) - |52) - 2|59) - 2|60) 

+2|61) + 2|62) - 8|79) + 8|82) - 4|84), 
H\7l) = 2|37) + 2|38) + 2|43) - 4|46) - 8|78) + 4|84), 
H\72) = -4|48) + |49) + |50) + |51) + |52) + 2|59) + 2|60) - 2|61) - 2|62) 

-8|81)-8|82), 
if |73) = -|49) - |50) - |51) - |52) - 4|54) + 2|59) + 2|60) - 2|61) - 2|62) 

+4|84)-8|85), 
iJ|74) = 2|35) + |50) - |52) + 2|55) + 2|56) - |57) - |58) - 2|61) + 2|62) + 4|83), 
H\75) = 2|35) + 2|42) - |50) + |52) + 2|55) - |57) + |58) - 2|61) + 2|62) + 4|83), 
H\76) = -|42) - |50) + |52) + |56) + |58), 
H\77) = |67), 

H\78) = |66) - |67) - |71), 
H\79) = |67) - |69) - |70), 
^|80) = -|68)- |69), 
H\81) = |68)-|72), 
H\82) = |70) - |72), 
^|83) = |74) + |75), 
^|84) = -|70) + |71) + |73), 
iJ|85) = -|66) - |69) - |73). 



(Al) 
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APPENDIX B: THE U INDEPENDENT EIGENVECTORS CORRESPONDING 
TO ZERO EIGENVALUE. 

The U independent eigenvalues corresponding to zero eigenvalues are the following 

^i) = |35) + |57), 

O2) = |55) + |57)-|65), 

^3) = |37)-|38) + |51)-|49), 

e^) = |60) + |61) - |63) + |64) - |65) + |83), 

65) = |60) + |62)-|63) + |64)-|83), 

Oq) = |77) + |78) - |80) - |81) + |82) + |84) + |85), 

67) = -|77) + |79) - |80) - |81) + |82), 

Os) = |10)-|11) + |49)-|52) + |58), 

^9) = |13) - |14) + |63) - |64) + |55) - |83), 

^10) = -|5) + |15) + |44) - |47) - |53) - |77) + |80), 

^11) = |14) - |15) + |46) - |48) - |60) + |84), 

612) = |44) + |45) - |47) - |48) + |53) - 2|77) + |79) + |81), 

9^3) = |29) - |28) + |30) - |31) + |75) - |74) + 2|33) - 2|34) - 2|76), 

Ou) = -|45) + |46) + |54) + 2|77) - 2|53) - 2|79) + 2|80) + |84), 

^15) = |37) - |44) + |47) - |49) - |77) + |53) + |80) - |82) - |84), 

^i6) = |43)-2|37) + 2|49)-|50)-|52), 

6^7) = |10) - 2|12) - |15) - 2|44) - |45) - |48) - |52) + |54) + |58) + 2|77) + 2|80), 

9^s) = -2|5) - 2|15) + 2|10) + 2|11) - 8|12) - |43) - 8|47) - 8|48) - 2|49) + |50) 

-|52) + 4|54) + 2|58) + 4|77) + 12|80) + 8|81) - 2|82) + 2|84), (Bl) 

For examplification we present also two U independent eigenfunctions corresponding to 
non-zero eigenvalues, namely +8 (the vector |fi)), and —8 (the vector |f2)). 

\vi) = |5) - |15) + |23) + |26) + |27) + |37) + |38) + |43) + 2|44) - |46) - 2|47) + |48) 
-|49) - |50) - |51) - |52) - 2|53) - |54) + 2|67) + 2|68) + 2|69) + |71) - |72) 
+ |73) + 4|77) - 2|78) - 4|80) + 2|81) + |82) + |84) - 2|85), 
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\V2) = |5) - |15) - |23) - |26) - |27) + |37) + |38) + |43) + 2|44) - |46) - 2|47) + |48) 
-|49) - |50) - |51) - |52) - 2|53) - |54) - 2|67) - 2|68) - 2|69) - |71) + |72) 
-|73) + 4|77) - 2|78) - 4|80) + 2|81) + |82) + |84) - 2|85), (B2) 
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FIG. 1: The site numbering for a 4 x 4 cluster, and the possible particle displacements for four 
electrons. 
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FIG. 2: The effect of the translation operator T. The result in the right hand side contains 16 
terms, from which, the first 8 are ploted in the figure. 



21 



T 



ra^0 







.0 









"^''- 

s, 

^ N 



+ 



FIG. 3: The effect of tlie translation operator T leading to sign changes in the mathematical 
expressions. 



22 



1) = r 



T 



|2> =T 



+ 



-T 



|3} =T 



+ 



+ 



+ 



-T 



+ 



+ 



|4) = \T 



2-^ 



|5} =T 



X 







\6)=T 



+ 



+ 



+ 



-T 











+ 








+ 








+ 
4 








\ 








1 








1 


1 










r 


K 


\ 


p 


r 




1 s 


1 


P 


► 







FIG. 4: Base wave vectors |1) 
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FIG. 5: Base wave vectors |7) — |11). 
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FIG. 6: Base wave vectors |12) — |18). 
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FIG. 7: Base wave vectors |19) - |23). 
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FIG. 8: Base wave vectors |24) - |27). 
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FIG. 10: Base wave vectors |32) - |38). 
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FIG. 12: Base wave vectors |46) - |50). 
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FIG. 13: Base wave vectors |51) — |56). 
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FIG. 14: Base wave vectors |57) - |62). 
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FIG. 15: Base wave vectors |63) — |69). 
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FIG. 16: Base wave vectors |70) — |73). 
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FIG. 17: Base wave vectors |74) - |79). 
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FIG. 18: Base wave vectors |80) - |85). 
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